3.777 \(\int \frac{x^3}{\left (a+b x^2\right )^2 \left (c+d x^2\right )^{5/2}} \, dx\)

Optimal. Leaf size=170 \[ \frac{a}{2 b \left (a+b x^2\right ) \left (c+d x^2\right )^{3/2} (b c-a d)}+\frac{3 a d+2 b c}{2 \sqrt{c+d x^2} (b c-a d)^3}+\frac{3 a d+2 b c}{6 b \left (c+d x^2\right )^{3/2} (b c-a d)^2}-\frac{\sqrt{b} (3 a d+2 b c) \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{c+d x^2}}{\sqrt{b c-a d}}\right )}{2 (b c-a d)^{7/2}} \]

[Out]

(2*b*c + 3*a*d)/(6*b*(b*c - a*d)^2*(c + d*x^2)^(3/2)) + a/(2*b*(b*c - a*d)*(a +
b*x^2)*(c + d*x^2)^(3/2)) + (2*b*c + 3*a*d)/(2*(b*c - a*d)^3*Sqrt[c + d*x^2]) -
(Sqrt[b]*(2*b*c + 3*a*d)*ArcTanh[(Sqrt[b]*Sqrt[c + d*x^2])/Sqrt[b*c - a*d]])/(2*
(b*c - a*d)^(7/2))

_______________________________________________________________________________________

Rubi [A]  time = 0.492164, antiderivative size = 170, normalized size of antiderivative = 1., number of steps used = 6, number of rules used = 5, integrand size = 24, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.208 \[ \frac{a}{2 b \left (a+b x^2\right ) \left (c+d x^2\right )^{3/2} (b c-a d)}+\frac{3 a d+2 b c}{2 \sqrt{c+d x^2} (b c-a d)^3}+\frac{3 a d+2 b c}{6 b \left (c+d x^2\right )^{3/2} (b c-a d)^2}-\frac{\sqrt{b} (3 a d+2 b c) \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{c+d x^2}}{\sqrt{b c-a d}}\right )}{2 (b c-a d)^{7/2}} \]

Antiderivative was successfully verified.

[In]  Int[x^3/((a + b*x^2)^2*(c + d*x^2)^(5/2)),x]

[Out]

(2*b*c + 3*a*d)/(6*b*(b*c - a*d)^2*(c + d*x^2)^(3/2)) + a/(2*b*(b*c - a*d)*(a +
b*x^2)*(c + d*x^2)^(3/2)) + (2*b*c + 3*a*d)/(2*(b*c - a*d)^3*Sqrt[c + d*x^2]) -
(Sqrt[b]*(2*b*c + 3*a*d)*ArcTanh[(Sqrt[b]*Sqrt[c + d*x^2])/Sqrt[b*c - a*d]])/(2*
(b*c - a*d)^(7/2))

_______________________________________________________________________________________

Rubi in Sympy [A]  time = 47.5891, size = 141, normalized size = 0.83 \[ - \frac{a}{2 b \left (a + b x^{2}\right ) \left (c + d x^{2}\right )^{\frac{3}{2}} \left (a d - b c\right )} - \frac{\sqrt{b} \left (\frac{3 a d}{2} + b c\right ) \operatorname{atan}{\left (\frac{\sqrt{b} \sqrt{c + d x^{2}}}{\sqrt{a d - b c}} \right )}}{\left (a d - b c\right )^{\frac{7}{2}}} - \frac{\frac{3 a d}{2} + b c}{\sqrt{c + d x^{2}} \left (a d - b c\right )^{3}} + \frac{\frac{3 a d}{2} + b c}{3 b \left (c + d x^{2}\right )^{\frac{3}{2}} \left (a d - b c\right )^{2}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  rubi_integrate(x**3/(b*x**2+a)**2/(d*x**2+c)**(5/2),x)

[Out]

-a/(2*b*(a + b*x**2)*(c + d*x**2)**(3/2)*(a*d - b*c)) - sqrt(b)*(3*a*d/2 + b*c)*
atan(sqrt(b)*sqrt(c + d*x**2)/sqrt(a*d - b*c))/(a*d - b*c)**(7/2) - (3*a*d/2 + b
*c)/(sqrt(c + d*x**2)*(a*d - b*c)**3) + (3*a*d/2 + b*c)/(3*b*(c + d*x**2)**(3/2)
*(a*d - b*c)**2)

_______________________________________________________________________________________

Mathematica [A]  time = 0.731788, size = 137, normalized size = 0.81 \[ \frac{1}{6} \left (\frac{\sqrt{c+d x^2} \left (\frac{6 (a d+b c)}{c+d x^2}+\frac{2 c (b c-a d)}{\left (c+d x^2\right )^2}+\frac{3 a b}{a+b x^2}\right )}{(b c-a d)^3}-\frac{3 \sqrt{b} (3 a d+2 b c) \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{c+d x^2}}{\sqrt{b c-a d}}\right )}{(b c-a d)^{7/2}}\right ) \]

Antiderivative was successfully verified.

[In]  Integrate[x^3/((a + b*x^2)^2*(c + d*x^2)^(5/2)),x]

[Out]

((Sqrt[c + d*x^2]*((3*a*b)/(a + b*x^2) + (2*c*(b*c - a*d))/(c + d*x^2)^2 + (6*(b
*c + a*d))/(c + d*x^2)))/(b*c - a*d)^3 - (3*Sqrt[b]*(2*b*c + 3*a*d)*ArcTanh[(Sqr
t[b]*Sqrt[c + d*x^2])/Sqrt[b*c - a*d]])/(b*c - a*d)^(7/2))/6

_______________________________________________________________________________________

Maple [B]  time = 0.023, size = 2400, normalized size = 14.1 \[ \text{result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  int(x^3/(b*x^2+a)^2/(d*x^2+c)^(5/2),x)

[Out]

5/12/b*a*d/(a*d-b*c)^2/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(1/2)/b*(x-1/b*(-a*b
)^(1/2))-(a*d-b*c)/b)^(3/2)+1/4/b^2*(-a*b)^(1/2)/(a*d-b*c)/(x-1/b*(-a*b)^(1/2))/
((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(
3/2)+1/2/(a*d-b*c)^2/((x+1/b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^
(1/2))-(a*d-b*c)/b)^(1/2)+1/2/(a*d-b*c)^2/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(
1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2)+1/2/b/(a*d-b*c)^2*(-a*b)^(1/2)/c/
((x+1/b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(
1/2)*x*d+1/2/b^2*d*(-a*b)^(1/2)/(a*d-b*c)/c/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)
^(1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(3/2)*x+1/b^2*d*(-a*b)^(1/2)/(a*d-b*c
)/c^2/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c
)/b)^(1/2)*x-1/2/b/(a*d-b*c)^2*(-a*b)^(1/2)/c/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*
b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2)*x*d-1/b^2*d*(-a*b)^(1/2)/(a*d
-b*c)/c^2/((x+1/b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d
-b*c)/b)^(1/2)*x-5/4*a*d/(a*d-b*c)^3/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(1/2)/
b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2)-5/4*a*d/(a*d-b*c)^3/((x+1/b*(-a*b)^(1/
2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2)-1/6/b/(a*d-b*
c)/((x+1/b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b
)^(3/2)-1/2/(a*d-b*c)^2/(-(a*d-b*c)/b)^(1/2)*ln((-2*(a*d-b*c)/b-2*d*(-a*b)^(1/2)
/b*(x+1/b*(-a*b)^(1/2))+2*(-(a*d-b*c)/b)^(1/2)*((x+1/b*(-a*b)^(1/2))^2*d-2*d*(-a
*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2))/(x+1/b*(-a*b)^(1/2)))-1/6/b
/(a*d-b*c)/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*
d-b*c)/b)^(3/2)-1/2/(a*d-b*c)^2/(-(a*d-b*c)/b)^(1/2)*ln((-2*(a*d-b*c)/b+2*d*(-a*
b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))+2*(-(a*d-b*c)/b)^(1/2)*((x-1/b*(-a*b)^(1/2))^2*d
+2*d*(-a*b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2))/(x-1/b*(-a*b)^(1/2)
))-5/12/b^2*(-a*b)^(1/2)*d^2*a/(a*d-b*c)^2/c/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b
)^(1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(3/2)*x+5/4*a*d/(a*d-b*c)^3/(-(a*d-b
*c)/b)^(1/2)*ln((-2*(a*d-b*c)/b+2*d*(-a*b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))+2*(-(a*d
-b*c)/b)^(1/2)*((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))
-(a*d-b*c)/b)^(1/2))/(x-1/b*(-a*b)^(1/2)))+5/12/b*a*d/(a*d-b*c)^2/((x+1/b*(-a*b)
^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(3/2)+5/4*a*d/(
a*d-b*c)^3/(-(a*d-b*c)/b)^(1/2)*ln((-2*(a*d-b*c)/b-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a
*b)^(1/2))+2*(-(a*d-b*c)/b)^(1/2)*((x+1/b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(
x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2))/(x+1/b*(-a*b)^(1/2)))-1/4/b^2*(-a*b)^(1/
2)/(a*d-b*c)/(x+1/b*(-a*b)^(1/2))/((x+1/b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(
x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(3/2)-1/2/b^2*d*(-a*b)^(1/2)/(a*d-b*c)/c/((x+1/
b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(3/2)*x
+5/4/b*(-a*b)^(1/2)*d^2*a/(a*d-b*c)^3/c/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(1/
2)/b*(x-1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2)*x-5/6/b^2*(-a*b)^(1/2)*d^2*a/(a*d-b
*c)^2/c^2/((x-1/b*(-a*b)^(1/2))^2*d+2*d*(-a*b)^(1/2)/b*(x-1/b*(-a*b)^(1/2))-(a*d
-b*c)/b)^(1/2)*x+5/6/b^2*(-a*b)^(1/2)*d^2*a/(a*d-b*c)^2/c^2/((x+1/b*(-a*b)^(1/2)
)^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(1/2)*x-5/4/b*(-a*b)^
(1/2)*d^2*a/(a*d-b*c)^3/c/((x+1/b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-
a*b)^(1/2))-(a*d-b*c)/b)^(1/2)*x+5/12/b^2*(-a*b)^(1/2)*d^2*a/(a*d-b*c)^2/c/((x+1
/b*(-a*b)^(1/2))^2*d-2*d*(-a*b)^(1/2)/b*(x+1/b*(-a*b)^(1/2))-(a*d-b*c)/b)^(3/2)*
x

_______________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \[ \text{Exception raised: ValueError} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x^3/((b*x^2 + a)^2*(d*x^2 + c)^(5/2)),x, algorithm="maxima")

[Out]

Exception raised: ValueError

_______________________________________________________________________________________

Fricas [A]  time = 0.405653, size = 1, normalized size = 0.01 \[ \text{result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x^3/((b*x^2 + a)^2*(d*x^2 + c)^(5/2)),x, algorithm="fricas")

[Out]

[-1/24*(3*((2*b^2*c*d^2 + 3*a*b*d^3)*x^6 + 2*a*b*c^3 + 3*a^2*c^2*d + (4*b^2*c^2*
d + 8*a*b*c*d^2 + 3*a^2*d^3)*x^4 + (2*b^2*c^3 + 7*a*b*c^2*d + 6*a^2*c*d^2)*x^2)*
sqrt(b/(b*c - a*d))*log((b^2*d^2*x^4 + 8*b^2*c^2 - 8*a*b*c*d + a^2*d^2 + 2*(4*b^
2*c*d - 3*a*b*d^2)*x^2 + 4*(2*b^2*c^2 - 3*a*b*c*d + a^2*d^2 + (b^2*c*d - a*b*d^2
)*x^2)*sqrt(d*x^2 + c)*sqrt(b/(b*c - a*d)))/(b^2*x^4 + 2*a*b*x^2 + a^2)) - 4*(3*
(2*b^2*c*d + 3*a*b*d^2)*x^4 + 11*a*b*c^2 + 4*a^2*c*d + 2*(4*b^2*c^2 + 8*a*b*c*d
+ 3*a^2*d^2)*x^2)*sqrt(d*x^2 + c))/(a*b^3*c^5 - 3*a^2*b^2*c^4*d + 3*a^3*b*c^3*d^
2 - a^4*c^2*d^3 + (b^4*c^3*d^2 - 3*a*b^3*c^2*d^3 + 3*a^2*b^2*c*d^4 - a^3*b*d^5)*
x^6 + (2*b^4*c^4*d - 5*a*b^3*c^3*d^2 + 3*a^2*b^2*c^2*d^3 + a^3*b*c*d^4 - a^4*d^5
)*x^4 + (b^4*c^5 - a*b^3*c^4*d - 3*a^2*b^2*c^3*d^2 + 5*a^3*b*c^2*d^3 - 2*a^4*c*d
^4)*x^2), 1/12*(3*((2*b^2*c*d^2 + 3*a*b*d^3)*x^6 + 2*a*b*c^3 + 3*a^2*c^2*d + (4*
b^2*c^2*d + 8*a*b*c*d^2 + 3*a^2*d^3)*x^4 + (2*b^2*c^3 + 7*a*b*c^2*d + 6*a^2*c*d^
2)*x^2)*sqrt(-b/(b*c - a*d))*arctan(-1/2*(b*d*x^2 + 2*b*c - a*d)/(sqrt(d*x^2 + c
)*(b*c - a*d)*sqrt(-b/(b*c - a*d)))) + 2*(3*(2*b^2*c*d + 3*a*b*d^2)*x^4 + 11*a*b
*c^2 + 4*a^2*c*d + 2*(4*b^2*c^2 + 8*a*b*c*d + 3*a^2*d^2)*x^2)*sqrt(d*x^2 + c))/(
a*b^3*c^5 - 3*a^2*b^2*c^4*d + 3*a^3*b*c^3*d^2 - a^4*c^2*d^3 + (b^4*c^3*d^2 - 3*a
*b^3*c^2*d^3 + 3*a^2*b^2*c*d^4 - a^3*b*d^5)*x^6 + (2*b^4*c^4*d - 5*a*b^3*c^3*d^2
 + 3*a^2*b^2*c^2*d^3 + a^3*b*c*d^4 - a^4*d^5)*x^4 + (b^4*c^5 - a*b^3*c^4*d - 3*a
^2*b^2*c^3*d^2 + 5*a^3*b*c^2*d^3 - 2*a^4*c*d^4)*x^2)]

_______________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \[ \text{Timed out} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x**3/(b*x**2+a)**2/(d*x**2+c)**(5/2),x)

[Out]

Timed out

_______________________________________________________________________________________

GIAC/XCAS [A]  time = 0.237834, size = 351, normalized size = 2.06 \[ \frac{\frac{3 \, \sqrt{d x^{2} + c} a b d^{2}}{{\left (b^{3} c^{3} - 3 \, a b^{2} c^{2} d + 3 \, a^{2} b c d^{2} - a^{3} d^{3}\right )}{\left ({\left (d x^{2} + c\right )} b - b c + a d\right )}} + \frac{3 \,{\left (2 \, b^{2} c d + 3 \, a b d^{2}\right )} \arctan \left (\frac{\sqrt{d x^{2} + c} b}{\sqrt{-b^{2} c + a b d}}\right )}{{\left (b^{3} c^{3} - 3 \, a b^{2} c^{2} d + 3 \, a^{2} b c d^{2} - a^{3} d^{3}\right )} \sqrt{-b^{2} c + a b d}} + \frac{2 \,{\left (3 \,{\left (d x^{2} + c\right )} b c d + b c^{2} d + 3 \,{\left (d x^{2} + c\right )} a d^{2} - a c d^{2}\right )}}{{\left (b^{3} c^{3} - 3 \, a b^{2} c^{2} d + 3 \, a^{2} b c d^{2} - a^{3} d^{3}\right )}{\left (d x^{2} + c\right )}^{\frac{3}{2}}}}{6 \, d} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x^3/((b*x^2 + a)^2*(d*x^2 + c)^(5/2)),x, algorithm="giac")

[Out]

1/6*(3*sqrt(d*x^2 + c)*a*b*d^2/((b^3*c^3 - 3*a*b^2*c^2*d + 3*a^2*b*c*d^2 - a^3*d
^3)*((d*x^2 + c)*b - b*c + a*d)) + 3*(2*b^2*c*d + 3*a*b*d^2)*arctan(sqrt(d*x^2 +
 c)*b/sqrt(-b^2*c + a*b*d))/((b^3*c^3 - 3*a*b^2*c^2*d + 3*a^2*b*c*d^2 - a^3*d^3)
*sqrt(-b^2*c + a*b*d)) + 2*(3*(d*x^2 + c)*b*c*d + b*c^2*d + 3*(d*x^2 + c)*a*d^2
- a*c*d^2)/((b^3*c^3 - 3*a*b^2*c^2*d + 3*a^2*b*c*d^2 - a^3*d^3)*(d*x^2 + c)^(3/2
)))/d